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Zusammenfassung 

Wir zeigen, daß die Grundzustandsenergie Em{Z) des Müllerfunktio- 
nals von (neutralen) Atomen der Ordnungszahl Z mit der quantenmecha- 
nischen Grundzustandsenergie Es{Z) bis zur Ordnung o{Z^''^) überein- 
stimmt, d.h. 

Em{Z) ^Es{Z)+o{Z3). 

1 Einleitung 

1.1 Quantenmechanik 

Der Zustandsraum eines iV-Elektronensystems sind die Einheitsstrahlen im iV- 
Teilchen- Hilbertraum 

N 
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Der Schrödingeroperator eines iV-Elektronenatoms 



= f-A„ - ^) + 5: (1) 

ri=l ^ II I / l<m<n<N ' ' 

ist in diesem Raum kanonisch als Störung des Laplaceoperators selbstadjungiert 
realisiert. Dessen Grunzustandsenergie ist Es{N,Z) := mi{a{SN,z))- Wir ver- 
merken für später, daß wir im neutralen Fall, also N = das erste Argument 
aller auftauchenden Energiefunktionale einfach unterdrücken, d.h. wir werden 
z.B. kurz Es{Z) für Es{Z,Z) schreiben. 



1.2 Funktional der Einteilchendichtematrix 

Wir nennen 7 eine (fermionische) Einteilchendichtematrix, wenn 7 ein selbstad- 
jungierter Operator der Spurklasse 6^(i^(K3) ® C^) auf ^^(R^) ® ist und 
< 7 < 1 gilt. Die Spur von 7 heißt Teilchenzahl. Wir bezeichnen die Menge 
der Einteilchendichtematrizes mit 

2::= {7 e 6i(i2(R3) ^C^) I < 7 < 1} (2) 

und setzen 

In -{lell Sp7<iV}. (3) 

Nach Gilbert I14j kann die Grundzustandsenergie eines jeden quantenmecha- 
nischen Coulombsystems - analog zum Hohenberg-Kohn-Theorem - als Mini- 
mum eines universellen Funktionais der Einteilchendichtematrix 7 geschrieben 
werden. 

Obwohl dieses Funktional unbekannt ist, gibt es mehrere Funktionale, die ei- 
ne gute Approximation für Es{Z) liefern. Das bekannteste ist wohl das Hartree- 
Fock- Funktional in der von Lieb |2_0j angegebenen Form. Um das Funktional ein- 
zuführen, verwenden wir folgende Notation: Zu gegebenem 7 G I sei (/^i, (^2, • ■ • 
ein vollständiges System orthonormierter Eigenvektoren mit zugehörigen Eigen- 
werten Ai, A2, . . . . Dann ist 7(x, y) := J2'^=o ^n^n{x)^n{y) der Integralkern von 
7 und p~f{i) :— Y^(j^i"f{x,x), die spinsummierte Diagonale, die Teilchendichte 
von 7. (Wie nicht unüblich, bezeichnen wir mit x und y Raum-Spin- Variablen, 
also x^y e r := R'^ x {1,2} und später mit Jp dx die Integration bezüglich j: 
und Summation über er.) Ferner sei 

das Coulombskalarprodukt und mit X das Austausch-Skalarprodukt 



2Jr Jt 

(Die zugehörigen quadratischen Formen bezeichnen wir mit denselben Buchsta- 
ben, setzen jedoch das Argument in eckige Klammern.) 
Damit können wir das Hartree-Fock-Funktional 

£nF-1 ^ K (4) 

7 ^ M{-^-z/\-M + D[p.^-x[^] (5) 
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definieren. Die Grundzustandsenergie des Hartree- Fock- Funktionais, kurz die 
Hartree-Fock-Energie, ist 

^HF(iV, Z) = inf{fHF(7) \i&^n}- (6) 

Es ist bekannt, dai3 das Infimum in ([ü]) für ganzzahliges N unverändert bleibt, 
wenn man die Minimierung auf idempotente Dichtematrizes, also solche, die von 
Slaterdeterminanten herrühren, beschränkt (Lieb [5Üj). Ferner ist bekannt, daß 
das Infimum für N < Z + \ stets für eine Dichtematrix mit Teilchenzahl N 
angenommen wird, d.h. ein Grundzustand mit Teilchenzahl N existiert (Lieb 
und Simon [HHg). 

Da für idempotentes 7 — + • • • + \lpis!){(Pn\ niit orthonormalen Spi- 

noren ipi, - ■ ■ , ipiq die Identität 

'^^hf(7) = N\^'^{ipi A ■ • • A ipN, SN.zfi A ■ • ■ A (Pn) 

gilt, ist die Hartree-Fock-Energie Es{N,Z) stets eine obere Schranke an die 
quantenmechanische Grundzustandsenergie. 

Ein weiteres Funktional der Einteilchendichtcmatrix wurde von Müller [25] 
angegeben 

f M : 2: ^ K (7) 
7 ^ Sp[{-A-Z/\-\h] + D[p^]-X[j^] (8) 
Die Müllerenergie ist 

EuiN, Z) = inf{fM(7) I 7 e Iat}. (9) 

Frank u.a. [T^ zeigten, dal3 für N < Z der - auch ohne diese Einschränkung 
existierende Minimierer - stets die Teilchenzahl N besitzt, d.h. es exisitiert 
ein Grundzustand mit Teilchenzahl N. Im Gegensatz zu einem Hartree-Fock- 
Minimicrer 7hf, dessen Bildraum stets N dimensional ist, ist der Bildraum ei- 
nes MüUerminimierers jedoch immer unendlich-dimcnsional. Offensichtlich ist, 
wiederum per constructionem, daß die Müllerenergie kleiner als die Hartree- 
Fock-Energie ist. Weiter gilt, daß im Zweielektronenfall (auch mehrere Zentren 
sind erlaubt), die Müllerenergie sogar die quantenmechanische Energie nach un- 
ten beschränkt [12], ein Faktum, daß durch numerische Rechnungen für alle 
Elektronenzahlen N nahe gelegt wird; mehr noch, diese numerischen Resultate 
deuten darauf hin, daß die Müllerenergien nicht mehr von den quantenmechani- 
schen Energien abweichen als die Hartree-Fock-Energien. Während jedoch auf 
der analytischen Seite für die Hartree-Fock-Energie 

EMZ)^EsiZ)+o{zl) (10) 

bekannt ist (Bach [2]) - dieses Resultat gilt auch für Moleküle und Ionen, wenn 
die Neutralität des Systems nicht zu stark verletzt ist -, fehlt ein entsprechen- 
des Resultat für die Müllerenergie. Es könnte also sein, daß die Müllerenergie 
unrealistisch klein ist. Daß dieses nicht zutrifft, ist der Inhalt der vorliegenden 
Arbeit. Genauer gesagt ist unser Resultat 

Satz 1. Es gibt eine positive Konstante c so, daß für alle Z > 1 
Em{Z) < EiiF{Z) < Em{Z) + czf- liö 

gilt. 
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Wir bemerken, daß dieses Resultat auf Moleküle mit nicht zu kleinen Kern- 
abständen - einschließlich der physikalischen Abstände - übertragbar ist. Ferner 
ist die strikte Neutralität nicht wesentlich. Das System muß nur im wesentlichen 
neutral sein. Wir präsentieren jedoch der Übersichtlichkeit halber hier nur den 
atomaren Fall. 

Weiterhin sei bemerkt, daß eine durch numerische Resultate nahe gelegte 
Ungleichung Es{Z) > Em{Z) auch sofort das Bachsche Resultat pU)) implizieren 
würde. 

Desweiteren möchten wir darauf verweisen, daß erst jüngst Dichtefunktionale 
hergeleitet und untersucht wurden, die die gleiche asymptotische Genauigkeit 
(bis zur Ordnung o{Z^^^)) besitzen, wie das Müllerfunktional (Lee u. a. [16]). 
Während die Konstruktion solcher Funktionale erheblichen Aufwand erfordert, 
ist dieses im Müllerfunktional ohne weiteres eingebaut. 

Weiter bemerken wir, daß Satz[T]gültig bleibt, wenn man im Müllerfunktional 
X[7^/^] durch J'ir[7"] mit 1/2 < a < 1 ersetzt. Solche Funktionale wurden von 
Sharma u.a. [27j eingeführt. 

Schließlich impliziert der Satz [T] mittels des obigen Bachschen Resultats 
und der asymptotischen Entwicklung der quantenmechanischen Grundszustand- 
senergie (Fefferman und Seco [ilinillllllllllSiaiH]) sofort, daß 

Em{Z) ^ Es{Z) + o{zi) = Etf{Z) + \z'^ - c<s,z\ + o(zt) (11) 

gilt. Hier ist E'^^iZ) die Thomas- Fermi- Energie des neutralen Atoms der Ord- 
nungszahl Z und CS ist die Schwingerkonstante (Schwinger |2Bj). 

Der Beweis des Satzes [T] ist von der Bachschen Idee [T] geprägt, die Ener- 
giedifferenz zwischen zwei Funktionalen durch eine Norm der abgeschnittenen 
Einteilchendichtematrix 7(1 — 7) abzuschätzen und diese Norm dann als klein 
zu erkennen. Dieses wiederum wird durch Vergleich mit einem geeigneten ap- 
proximativen Zustand, hier dem halbklassichen Grundzustand des Schrödinge- 
roperators mit Thomas-Fermi-Potential, gewonnen, wobei unser Vorgehen von 
Graf und Solovej [15j inspiriert ist. 

2 Einfache Schranke an die Müllerenergie 

5 

Wir benötigen eine einfache Schranke der Größenordnung Z'i an die Mülle- 
renergie, um die Asymptotik in dieser Größenordnung zu finden. Eine solche 
Schranke wird in diesem Abschnitt entwickelt. 

Cioslowski und Pernal [3] (siehe auch [T^) konnten zeigen, daß der MüUer- 
sche Austausch-Korrelations- Term keine einfache Ungleichung vom Lieb-Oxford- 

1 4 

Typ (Lieb [18], Lieb und Oxford [21]) der Form ^[72] < G / zuläßt Q Den- 
noch bewirkt dieser Term - wie der Austauschterm des Hartree-Fock-Funktio- 
nals - eine Energieerniedrigung, die höchstens der Größenordnung Z'i ist, was 
aus folgendem Lemma folgt; denn das reduzierte Hartree- Fock- Funktional 

frHF(7) Sp[(-A - Zl\ ■ 1)7] + D{p^\ (12) 

ist offensichtlich größer als das Hartree-Fock-Funktional, da ein negativer Term 
entfallen ist. Ferner sei i?rHF(A^, Z^ das Infimum von i?iHF auf Z^r. 

^Hier und im folgenden bezeichne C eine generischc Konstante. 
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Lemma 1. Sei 7 ein Minimierer des Müllerfunktionals auf In- Dann gilt 

Em{Z) < Ehf{Z) < E,hf{Z) < SrüFh) < Em{Z) + Czt. (13) 

Beweis. Die ersten drei Ungleichungen sind per constructionem und Variations- 
prinzip offensichtlich richtig, sodaß wir lediglich die letzte Schranke zeigen. Dazu 
wählen wir ein später näher zu bestimmendes e G (0,1) und bezeichnen mit 7 
einen Minimierer des Müllerfunktionals. Dann gilt 

ErMZ; Z) > EMZ, Z) > Em{Z., Z) - Sp[(-A - Z/\ ■ 1)7] + D[p^] - X[j^] 
> Sp[(-(1 - e) A - Z/1 • 1)7] + D[p,] + e Sp{-A^) - X[^^] 



> 



1 

^^E,MZ) + { ^ dx^ e|V,75(x,y)p - ^g^l^dy [7 £ In} 



>- 157 Jr I72 (a:,j/)Pdy 

> E,-aF{Z) - eCzl - e-^N/16 = E,-aF{Z) - Czt (14) 

Hier haben wir neben der Skalierung des reduzierten Hartree-Fock-Funktionals 

1 1 

benutzt, daß £'rHF(^) — E^f{^)Z^ + o(Z3 ) gilt, was aus dem Beweis von [TÜl 

Satz 5.1] folgt. Sodann beachten wir, daß N ~ Z gilt und wir deshalb e ~ Z 3 
wählen. □ 

Lemma 2 (Virialsatz (Frank u.a 12])). Sei 7 ein Minimierer des Müllerfunk- 
tionals £m auf In- Dann gilt Sp(— A7) = — £'m(7)- 

Beweis- Wir strecken den Grundzustand 7 und definieren für positives A die 
gestreckte Dichtematrix 7a durch ^\{x,y) :— \^j{X]:,a;\t),T). Die gestreckte 
Dichtematrix gehört offensichtlich auch zuXjv. Da /(A) := £m{i\) für positives 
Argument differenzierbar ist und für A = 1 ein Minimum hat, folgt die Aussage 
aus /'(l) = 0. □ 

7 7 

Da Efif{Z) = Etf{1)Z^ +o(Z3) (Lieb und Simon [24 ) gilt, folgt aus 
Lemmata [1] und [2 daß 

Sp(-A7) = 0{zl) (15) 

gilt. 



3 Verdampfungsgrad des Müllergrundzustandes 

Der Bildraum einer das Müllerfunktional minimierenden Eintcilchcndichtema- 
trix 7 ist unendlichdimensional (Frank u.a. |12|). eine Eigenschaft, die das 
Müllerfunktional mit der Schrödingertheorie teil (Friesecke jl3 , Lewin ^17]). 
Frank u.a. vermuten sogar, daß ein solches 7 trivialen Kern hat. Von der Ein- 
teilchendichtematrix eines quantenmechanischen Grundzustandes ist seit Bach 
[T] bekannt, daß sie nicht zu stark von einer Projektion abweichen kann. Ein 
Maß dafür ist die Spur der Dichtematrix 7(1 — 7), die verschwindet, wenn 7 
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eine Projektion ist. Eine stärkere Eigenschaft ist die fast vollständige Konden- 
sation der Dichtematrix in einen vorgegebenen Zustand, eine Eigenschaft, die im 
quantenmechanischen Fall von Graf und Solovej |.15j untersucht wurde. Prakti- 
kabilitäts halber wählen wir die halbklassiche Projektion P auf die gebundenen 
Zustände im Thomas-Fermi-Potential: Sei :— Z/\ • | — p * | • das Thomas- 
Fermi-Potential, das durch den Minimierer p des Thomas-Fermi-Funktionals 
Etp eines Atoms mit Ladungszahl Z gegeben ist. Ferner sei 



f(2^i?)-i|y|-isin(^|y|/i?) |y| < i? 
lo \i\>R 



mit R — Z ■^Z^. Weiter führen wir den kohärenten Zustand mit Implus p. Ort 
q und Spin er ein, d.h. 

|a) = |p,q,T) e'P-5(- - q)5.,, (16) 
mit a := (p, q, er) £ F := x x {1, 2}, und 



dn{a) := (27r)-3 / dp / dq ^ 

JE3 Jr3 2 



sei das natürliche Maß auf F im Sinne von Planck, das die Elektroncnzustände 
im Phasenraum zählt. Ferner sei M der klassisch erlaubte Bereich des Phasen- 
raumes, d.h. 

M:={(p,q,a) eF I p^ - 0(q) < 0}. 

Der „halbklassische Projektor"@ des Schrödingeroperators mit Thomas-Fermi- 
Potential ist 

dn{a)\a){a\. (17) 



M 



Damit haben wir 



Lemma 3. Sei 7 ein Grundzustand des Müllerfunktionals eines neutralen Atoms 
der Ordnungszahl Z . Dann gilt für den Verdampfungsgrad d{'j, P) des Zustands 
7 aus dem fermionischen Zustand P 

69 

(5(7,P) :=Sp((l-P)7) = 0(Z7Ö). (18) 
Beweis. Zunächst erinnern wir daran, daß für alle 7 G I 

Sp(7(^A^)) > ^ dr!(a)(p2 - 0(y))_ - CZ^-W (19) 

— ://tf 

gih (Lieb [Abschnitt V.A.2][IÜ], Thirring 28J). Damit haben wir 

frHF(7) > Sp(7i?TF) - D[p] > ^dr!(a)(p2 - 0(j:))_ - D[p] - CZ^^ , (20) 

da D[p — p-y,hf] — 0. Desweiteren benutzen wir die Phasenraumabschätzung 
(Bach [21 Lemma 12], Graf und Solovej Lemma 7 und den Beweis von 
(1.16)]) 

df^(a)(p2-0(j;)+a)=O(ai). (21) 

a<p2-0(j)<O 



^Es sei vermerkt, daß der so eingeführte Operator P kein Projektor ist. 
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Der Vollständigkeit halber fügen wir den Beweis der Abschätzung (|2ip von Graf 
und Solovej in leicht vereinfachter Form hier ein: 

/ dO(a)(p2 -(/.(y) +a) / db f dn{a) 

J-a<p2-0(f)<O Jo J-6<p2-</.(j)<0 

db / dy(0(y)2 - ~b)l)^2n db de (0(y) - c)l 



3 



6 /. 1 pa pb 

1-4 _ _ 
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<2n db de dy(C|yr^ - c) 2 = C / dfe / dce-I^Cei, (22) 
Jo ^0 Je-' Jo Jo 

wo wir das Potential (j) durch die Sommerfeldlösung C|j:|^'' nach oben abgeschätzt 
haben (Lieb O S. 607]). 
Weiter sei 

SmAi) ■= ^m{i) - aSp(l - P)7. 

Nun gilt 

aSp(l-P)7 = ^M(^)-fM.a(7)^-ErHF(^)-[frHF(7)-aSp((l-P)7)]+Czt 

(13) 

< Emf{Z) + aZ- Sp((i7TF + aP)7) + D[p] + Czt 
<^dr!(a)(p2-(^(y)+a)_-Z?[p]-^dr!(a)(p2-^(j:)+a;^A,(y))_+ö[p] + Czi 

69 7 69 

dr2(a)(p2 - 0(y) + a) + 0(Z3ö) < C(a4 + Zsö), (23) 



-a<p2-0(f)<O 
138 

d.h. für a = Z W5 erhalten wir 

69 

5(7,P)=0(Z7Ö), (24) 

also die behauptete Abschätzung an den Verdampfungsgrad. □ 

Wir möchten an dieser Stelle keine Vermutung über den tatsächlichen Ver- 
dampfungsgrad abgeben. Es ist aber klar, daß die hergeleitete Fehlerabschätzung 
nicht optimal ist. Benutzt man z.B. den Z^- Teilchen-Grundzustand von H^f und 
die von Fefferman und Seco |TÜ1 Hl [Hl [3 O El [S] angegebene asymptotische 

5 

Entwicklung, so ergibt sich ein Fehler der maximalen Größe 0{Z7). Wir hal- 
ten also fest, daß (jlSp zwar nicht optimal aber einfach und für unsere Zwecke 
hinreichend ist. 

Folgerung 1. Sei 7 ein Minimierer des Müllerfunktionals auf Tz- Dann ist 

Sp(7(l -7)) < Czfi. (25) 
Beweis. Da < 7 < 1 und Sp7 = SpP = Z ist, haben wir 

Sp(7(l-7)) = Sp(P7(l-7))+Sp((l-P)7(l-7)) < Sp(P(l-7))+Sp((l-P)7) 

69 

= Z - Sp(P7) + Sp((l - P)7) = 2 Sp((l - P)7) = 2<5(7, P) = 0{Z7Ö). (26) 

□ 
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4 Schranke an die Müllerenergie: Abgeschnitte- 
ne Einteilchendichtematrix 

Lemma 4. Sei 7 ein Minimierer des Müllerfunktionals Em auf 2z- Dann gilt 

< EnF{Z) ~ Em{Z) < Czi(Sp[7(l - 7)])5. 
Beweis. Gemäß Variationsprinzip haben wir 

< EMZ) ~EuiZ)<Uä.[ d,bi:M|^^M£2^ (27) 

2 ./r ./r r - tl 



1 1 
2./r Ir ^ \x-t)\ ^ ' 

1 



< 



dx I dy VSp|72(l -72)r (29) 



2^ 

< 2v'Sp(-AiVSp|7(l - 7)1 < Czl ^/Sp|7(l - t)1, (30) 

WO wir die Schwarzsehe Ungleiehung von ((28|) nach ([29|) und die Hardysche Un- 

1 

gleichung und 72 > 7 von P5)) nach ((^Ü)) sowie in der letzten Zeile Sp(— A7) = 
7 

0{Z3) (Gleichung ^) benutzt haben. □ 

Beweis des Satzes{^ Die untere Schranke ist trivial; denn die Hartree- Fock- 
Energie ist stets größer als die MüUerencrgie. Um die obere Schranke zu zeigen, 
müssen wir nun lediglich noch das Lemma [3] und die Folgerung [T] kombinieren 
und erhalten die behauptete Ungleichung des[TJ Satzes. □ 
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